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On the Theory of Linear Difference Equations.* 

By R. D. Cabmichael. 



Under the above general heading I present here a discussion of each of 
the following topics : 

I. On the Finite Integral of a Function. 

II. Entire Functions Defined by Linear Homogeneous Equations with 
Rational Coefficients. 
III. Applications of the Preceding Results to Non-homogeneous Equations. 

The three parts of the paper are closely connected by having the common 
purpose of investigating the character of functions defined by linear difference 
equations, especial emphasis being put on the possibility of obtaining funda- 
mental systems of solutions each function of which is an entire function. 
Where such a fundamental system does not exist, solutions are obtained having 
as simple sets of singularities as possible. 

Each part of the paper has its own introduction, containing references to 
the pertinent literature and giving a brief statement of the results obtained in 
that part. 

I. On the Finite Integbal of a Function". 

Introduction. 

Guichardf has proved that if G{x) is any entire function there exists 
another entire function u(x) satisfying the relation 

u(x + 1) — u{x) = G(x). 

Such a function u(x) he obtains in the form of an integral. AppellJ and 
Hurwitz § have given other proofs of the same theorem. Hurwitz, || by an 
elegant method, has proved the additional theorem that if <p(x) is any function 
of x which is meromorphic throughout the finite plane there exists a function 

* Read before the American Mathematical Society, Chicago, April 1912. 

f Annales Soientifiques de I'Eoole Normale Superieure, (3) 4 (1887) , pp. 361-380. 

J (Liouville) Journal de Hathematiques (4) 7 (1891), pp. 157-219; especially pp. 159-176. 

§Acta Mathematica 20 (1897), pp. 285-312. 

J| Loo. cit., pp. 308-312. 
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F(x) which is also meromorphic throughout the finite plane and satisfies the 

relation 

F{x + 1) — F(x) = Q(x). 

The primary object of the present part of this paper is to prove Guichard's 
important theorem by direct means. Accordingly, in §1 a power series for 
u{x) is assumed and direct substitution is made in the first equation above. 
There results an infinite system of equations for the determination of the 
coefficients in the series. This system is solved by what is believed to be a 
method of considerable value,* the solution obtained being such that the power 
series converges throughout the finite plane. 

In § 2 I determine properties of the finite integral of a function which is 
analytic outside of a given circle, except possibly at infinity. 

In §§3 and 4 I treat the same problem for functions which are analytic 
within a circle and within a circular ring respectively. 

§ 1. The Finite Integral of mi Entire Function. 
Let G(x) be any entire function and let it be expanded in the form 

G(x) = a Q -i- a 1 x + |y# 2 + -^jX 3 + 

We shall show that there exists an entire function u(x) satisfying the relation 

M(a> + 1) — u(x) = G(x). (1) 

Let us write 

u(x) = u + UjX + -|y ^ 2 + ^r xS + » 

and substitute in (1). If we equate the coefficients of like powers of x we 
have the following infinite system of equations in u lf w 2 , . . . ., 

2 -jy u i+i , = a t , * = 0,1,2, (2) 

It is clear that u is arbitrary. 

It is natural to expect that a solution of system (2) is linear in the a's. 
Accordingly, let us seek a solution in the form 

00 

u k =2a u a lt k = 1,2, (3) 

1 = 

If we substitute the value of u k from (3) into (2) we have 

00 -t 00 

X -jj Xa i+hl a l = a i , i = 0, 1, 2, 

j = l ] i 1 = 

* It is my purpose in a future paper to develop the general consequences of this method of dealing 
with non-homogeneous infinite systems of equations. 
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We shall determine the a's so that each of these equations is a formal identity 
in the a's. A necessary and sufficient condition for this is that the a's satisfy 
the doubly infinite system of equations 

.2 -jja i+i , , = S u , \l - ; x ; 2 ; ; ; ; ; ; (4) 

where S u is equal to 1 or according as i is or is not equal to I. 

For any fixed value of I we have the singly infinite system of equations 

X T ra i+j>l = S u , * = 0, 1, 2, (5) 

for determining those a's which have the given second subscript I. The 
(I + l)th equation of the set has the second member 1 ; all the other equations 
of the set have the second member 0. 

"We shall now obtain a particular solution a kl = a hl of (5) in which a M = 
when k > I -f- 1. For this case (5) reduces to the finite system (the order of 
equations being inverted) 

a i+i,i — 1> 
«*,*+ 2iai+ 1 ,i = 0| 

«i-i, i+ jj"i,i+ Jf^i+i, t = °> ( 6 ) 



«u + 2T«2i + + (T+T) | a i+i,i — °- 



This system may readily be solved by determinants, but a solution more con- 
venient for our purposes may be obtained in the form of a contour integral. 
If we write 

-L = c o + cj + c 2 P + c z e + • • . ., (7) 



e*- 



expand the denominator in powers of £, clear of fractions and in the resulting 
equation equate to zero the coefficients of the different powers of £, we have 
for determining the c's the following equations : 











Co 


= 1, 






c l 


+ 


Co 
2! 


= 0, 


c 2 


+ 


2! 


+ 


c 
3! 


= 0, 



166 Carmichael, : On the Theory of Linear Difference Equations. 



Comparing this system with (6) we see that a J+1 _^ ? = c A .. But from (7), by- 
means of the Cauchy theorem on residues, we have 

where the contour of integration G incloses the point zero and no other singu- 
larity of the integrand. Therefore 

1 r d£ 

«i+i-*,i = 2nV zZi Jc (e«— l)g* ' * = °i 1 » ••••i ? J a *i = °» £>* + !• 

Clearly this result may be written in the simpler form 

*"=d^rX (e*-iV-" * =1 ' 2,3 (8) 

If in (8) we let I run over all the numbers 0, 1, 2, . . . ., we obtain a com- 
plete solution of (4). Substituting from (8) into (3) we have 



Ui 



1 " c d % 

~2nV=l.=.°*Jc(e»-l)^+ 1 -*- 



(9) 



This result affords a formal solution of system (2). That in general it is 
divergent and hence does not afford an actual solution is easily shown. For 
we have 

£ _ i ;* c ** 



l.=o C Jo(e* — 1] 



e f -l 27tV^l*=<T Jo(e*-l)^' 

as one sees readily from the value of c k obtained above. Now the radius of 
convergence of this series is not greater than 2n, since £/ (e 1 — 1) has a pole at 
£ = 27i V — 1. Hence if a 4 = m\ as in the case when G (x) = e mx , where m>2n, 
series (9) diverges at least for Jc = 1. 

This (divergent) formal solution of (2) has been obtained by taking a 
particular solution of (4). It will now be shown that (4) has an infinitude of 
solutions, and that if an appropriately chosen particular solution be selected 
from this totality it will lead to an actual solution of (2). 

If we write a kl = a kl + {3 kl , substitute in (4) and make use of the fact that 
d kl is a solution of (4), it is clear that we have 

,1,7^".,.,= o, {?=8: i; 1: : :: ;: <i»> 

Thus for each value of I we have to solve a system of equations of the form 

2 -i- &+, = (), » = 0,1,2, (11) 



1=1 J' » 
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If ($(x) is any periodic function of period 1 which is analytic at zero, and 
we write 

£(20 = /3 + /^ + -ff-^+f!-^ + ----» 

substitute in the equation B(x + 1) — {3(x) = and equate to zero the coeffi- 
cients of like powers of x we have system (11). Therefore, for every periodic 
function of period 1, and analytic at zero, there exists a corresponding formal 
solution of (11) ; and this solution is furnished by the coefficients fi lf (3 2 , .... 
in the expansion of fi(x) above. If this expansion converges for x = 1 the 
formal solution becomes an actual solution. This result will be used in con- 
structing an appropriate solution of (10). 

From the function ce 2n " v ~ lx , where c is any constant and n is any positive 
or negative integer, we have a solution B k = c(2nnV — 1)* of (11). Now the 

residue of the function 

1 



at the point 2nnV— 1, is (2nnV— l)*-'- 1 . Hence, if c = (2nnV— 1) - 1 - 1 , we 
have @ k equal to the residue in consideration. That is, we have 



Pk 2nV=lJ (e*-l)£ | -*+ 1 ' 



where the contour of integration incloses the point 2nnV — 1 and no other 
singularity of the integrand. 

Now the sum of any finite number of solutions of system (11) is also a 
solution; hence it is clear that a set of values satisfying (10) is 

o 1 v C ^R [fc = 1, 2, . . . ., 

Pkl ~2nV=l J (e^-l^-H-!' 0=0,1,2,...., 

where for a fixed I the sum of 21 integrals is to be taken, one around each 
of the 21 contours inclosing separately each of the 21 points 2n7tV — 1, 
n = ± 1, ± 2, . . . . , ± I and no other singularity of the integrand. Combining 
this value of (3 U with that of a kl in (8) we have for a kl the value 



a 



kl 



2nV 



1 r d% p = l, 2, 

/ZTT J Cl (e«- l)|'-*+ 1 ' \l = 0, 1, 2, . . 



where the contour of integration C t incloses the points 2w7tV — 1, n = 0, 
±1, ±2, ...., ±1 and no other singularity of the integrand. 
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Substituting this value of a kl in (3) we have 
1 °° r d£ 

*°2W=5£ a 'J«. <*-l)V-+" k=1 ' 2 ' <12> 

In order to show that this is a valid solution of (2) we must first prove that 
the series is convergent for every value of Tc. 
Denote by l sh the integral 

Z «* = } Cs ( e «_i)£.-*+i ' (13 ) 

Choose as the contour of integration the square formed by four lines each at a 
distance (2s + l)n from the origin, two of them being parallel to the axis of 
reals and the other two being parallel to the axis of imaginaries. Making the 
transformation £ = (2s + l)nw, we have 



L,= 



1 r dw 

' sfc ~ (2« + l) s -*7i*- fc J s (e (2Hl)7 ™— l)w*~ k+1 ' 



where the contour 8 is the square formed by four lines each at a distance 

unity from the origin, two of them being parallel to the axis of reals and the 

other two being parallel to the axis of imaginaries. Along the contour of 

integration 1/ \e (2s+1 ^ nw — 1| remains finite; let us say that it is always less 

than A. Along the path of integration \w\ is equal to or greater than 1. 

Hence 

A c SA 

I 7 **' K (2s + l) s -*7i s - fc JJ^' = (2s + 1) '-*«*-* ' (14) 

From this result it follows that the series in (12) is convergent for every 
k provided that the same is true of the series 

QO 

V a > 

s =o(2s+l)*-*7l s - fc ' 
That this series is convergent follows at once by comparison with the con- 
vergent series «. 1 

stosl n°- k ' 

the convergence of the latter being a consequence of the fact that G (x) is 
entire. 

Substituting these values for u k in the expansion of u(x) and choosing a 
convenient value for the arbitrary constant u , we have the following result: 

wW =^bii.fTi„ a -X, («t_ &■-'+■ • <15) 

We shall next determine the circle of convergence of this power series. 
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Let R denote any circle of radius R about the point zero. We shall show 
that (15) is convergent in this circle R. Form the set of functions 

1 °° x k r d£ 
u ^ x) = YnV=ulokl a * X. R--ifF-^' S = °' lf 2 ' (16) 

In view of the inequality (14) it is easy to see that each of these series is 

convergent in the circle R. One has only to compare the 5th series with the 

series 

" R h 8A 

2 tt« 



*=oft!"" J (2s + 1) •-**•-* ' 

the convergence of the latter being obvious. 
Next we shall show that the series 

u (x) + Mi(a;) + w 2 (x)-j-... . (17) 

is uniformly convergent throughout R. By a theorem of Weierstrass it will 
then follow that the sum is analytic in R and may be written in the form of the 
second member of (15). From (14) and (16) we have 

\u,(x)\<a„ \x\<R, 

where 

1 " R h , , 8^ 



2n *=<,&! (2s + l)'~ k 7t°- k ' 

One shows readily that 

» Rk 

Ofcbt^la.-nl (2s + 1)' 1 »S,*T (2S + 3)kn * 
a. -~K| (2. + S)-P» jg (&+1)S(i 

- K+il (2g + i)* _! - 



Hence 



a. 



-.+i ^ |g.+il i J_ **. ^ Kti I i «"■" 



^ I *~«-H I _ g2«T < ^- 



\a s \ 2s+ 3 7t ^ \a,\ s + 1 2n ' 
Since G(x) is an entire function the series 

i\a,\±-R> 

s=0 S I 

is convergent for every value of R. Hence the ratio of two consecutive terms, 

i a i±i I X 5 
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is ultimately less than or equal to one. If R > e 2R7r /2n this ratio is greater 
than the ratio cc s+1 /a s ; and therefore the series 

a Q + a x + a 2 + a 3 + . . . . 
is convergent, since a s+1 /a s ultimately becomes and remains less than 1. But 
\ U A X )\ < a s> 1^1 = ^; an< ^ therefore series (16) is uniformly convergent 
throughout R. Hence its sum is analytic throughout R and can be written in 
the form (15). But R is any circle about the origin. Hence (15) converges 
throughout the entire plane; that is, u(x) is an entire function. 

It is easy to show directly that u(x) as defined in (15) actually satisfies 
relation (1). We have thus proved the following theorem: 

I. If G(x) is an entire function expansible in the form 

G(x) —a + a x x + -^y x 2 + -^- x 3 + , 

then there exist entire functions u(x) satisfying the difference equation 

u(x +1) — u(x) = G(x). 
One such entire function is 

1 * x k Z. r d% 



X) — ; Z ^-T 2, a s I 7-j — T r 



3-k + l ' 



where the contour of integration C s incloses the points 2nnV — 1, n = 0, ± 1, 
± 2, . . . ., ± s, and no other singularity of the integrand. 

Evidently, the general solution of the difference equation is obtained by 
adding to the above solution any arbitrary periodic function of period 1. If 
the solution is restricted to be entire this periodic function must be entire. 

The value of the integral I sJc , defined by equation (13), takes a very simple 
form when s — fc is even. For, then, the residues of the integrand at the two 
points 2nnV — 1 and — 2nnV — l,n=f=0, being respectively (2nnV — i)*-«-i 
and ( — 2nnV — l) k ~ s ~ 1 , are such that the part of the integral which comes 
from these two residues is zero. Hence 



= 1 



d£ 
.„ (e«-l)^-*+ 1 ' 



when s — & is even. 



That is, l sk is in this case 2tiV — 1 times the coefficient of £*-*+ 1 in the develop- 
ment of £/(e* — 1) in powers of £. Hence, when s — k = we have 
I sk = — 7i V — 1. It is easy to see that I sk is equal to zero when s — A; is even 
and not zero. To prove this it is sufficient to notice that 

-A-+U or i^±ll 
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is an even function; for then the coefficients of odd powers of £ (greater than 
the first) in the development of £/ ' (e i — 1) are all zero. 

Making use of these results as to the value of I sk we have at once the fol- 
lowing theorem : 

II. The function u(x) of Theorem I may be written in the form 

ivhere the accent on the second summation sign indicates that for each value 
of k only those values of s are to be taken for which s — k is odd. 

Corollary 1. If G(x) is an even function (so that a s is zero when s is 
odd), then the equation 

u(x + 1) — u{x) = G(x) 

has a solution which is the sum of — ^G(x) and an odd function which is 
entire. 

Corollary 2. // G{%) is an odd function, (so that a s is zero when s is 
even), then the equation 

u(x + 1) — u(x) = G(x) 

has a solution which is the sum of — ^G(x) and an even function which is 
entire. 

§ 2. The Finite Integral of a Function which is Analytic Outside 

of a Circle, except at Infinity. 

Let us consider the equation 

u(x) = G{x) + -?—+1>(x) = X(x), (18) 

Jb CI 

where G(x) is an entire function and $>{x) is analytic outside of a circle of 
finite radius R about the point a as a center, the development of <p(x) in nega- 
tive powers of x — a containing no constant term and no term in (x — a)~\ 
It is evident that a solution u(x) of (18) may be written in the form 

u(x) = u t (x) + n 2 (x) + u 3 (x), 

where u^x), u 2 (x), u 3 (x) are solutions respectively of the equations 

u 1 (x + 1) — u 1 (x) = G(x), 

(19) 



22 



u 2 (x-\- 1)- 


-uJx) = , 

2V ' x — a' 


u 3 (x+ 1)- 


-u 3 (x) = $(x) . 
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A solution u^x) of the first equation (19), analytic throughout the finite 
plane, was found in the preceding section. For the second and third equations 
we have the following two formal solutions : 

u i (x)= — + i(i- ^~W)> u 2 (x)=i( 1 , + i), 

x — a t=A * x — a -\- 1/ ' zv ' i=i\x — a — i t/' 

00 00 

u 3 (x)=—'Zq>(x + i), u z (x)—l l p{x — i). 

t=0 i=l 

Let us consider the first pair of formal solutions for u 2 (x) and u 3 (x). It is 
easy to see that there exists a line parallel to the axis of imaginaries to the 
right of which each of the series is uniformly convergent. Hence to the right 
of this line u 2 (x) and u 3 (x), as defined by these series, are analytic except at 
infinity. From this result and the fact that u x {x) is analytic throughout the 
finite plane it follows that u(x) is analytic to the right of the same line. Now, 
by means of the equation 

u(x — 1) = u(x) — X(x — 1), 

the value of u(x) is known at x — 1 when it is known at x . In this way the 
solution obtained above may be defined throughout any finite region of the 
plane. Thus we see that u(x) is analytic throughout the finite plane except at 
the singularities of X(x) in the finite plane and the points congruent to them 
on the left. (A point a — n is said to be congruent to a on the left if n is a 
positive integer ; • similarly, we say that a -\- n is congruent to a on the right 
if n is a positive integer.) 

In a similar way we may treat the second pair of formal solutions for 
u 2 (x) and u z {x). Here we consider first the part of the plane which is to the 
left of a suitably chosen line parallel to the axis of imaginaries. It is clear 
that the solution u(x) of (18), which we obtain in this case, is analytic 
throughout the finite plane except at the points congruent on the right to the 
singularities of X(x) in the finite plane. 

Thus we have the following theorem : 

III. Let X(x) be any function of x which is analytic outside of a circle 
about the center a of finite radius R, except possibly at infinity. Then there 
exists a function u(x) satisfying the relation 

u(x + 1) — u(x) = X(x) 

and having either one of the following two properties: 

1. It is analytic throughout the fintte plane except at the singularities of 
X(x) in the finite plane and the points congruent to them on the left. 

2. It is analytic throughout the finite plane except at points congruent 
on the right to the singularities of X(x) in the finite plane. 
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§3. The Finite Integral of a Function which is Analytic within a Circle. 

Let ^(x) be any function which is analytic within a circle R of radius R 
about the point a as a center. A method due to Weber * gives us in this case 
an immediate solution of the problem of finding a function u(x) which satisfies 
the equation 

u(x + 1) — u(x') = -^(tc) 

and is analytic in any circle within the circle JR. Let R be any circle about a 
as a center and having a radius less than that of the circle R. Denote by 
(RR) the circular ring between these two circles. On a line through a and 
parallel to the axis of imaginaries take two points A and B, lying within the 
circular ring (RR) , the imaginary coefficient of A being less than that of a and 
the imaginary coefficient of B being greater than that of a. Let x be any 
point in the circle R and' denote by I the line through the points . . . .x — 1, 

x, x -\- 1, 

Now form the integral 

«<->=X i- ( 2£-» i=v=l - < 20 > 

where the path of integration C extends from A to B, lies entirely within the 
circle B, passes through a point P lying on the line I between x — 1 and x and 
has no other point in common with the line I. Then, on account of the peri- 
odicity of the integrand considered as a function of x, it is clear that 

^(z)ds 



»(- + D = Xt3^' 



where the path of integration is as in the preceding case except that the point 
P is now between x and x + 1. The difference u(x + 1) — u(x) is expressible 
as a contour integral taken along the curve made up of C taken in a positive 
sense and G taken in a negative sense, the integrand being the same as above. 
Hence, by Cauchy's theorem, u(x-\-l) — u(x) = ^(x). But u(x), as defined 
by (20), is evidently an analytic function of x at every point x in the circle B. 
Hence we have the following theorem : 

IV. If "^{x) is analytic within a circle of radius B about the point a and 
R is any circle about a with radius less than B, then there exists a function 
u(x) which is analytic throughout R and satisfies the relation 

u(x + 1) — u(x) = ^(x). 

Such a function is defined in equation (20). 

*Acta Mathematica, Vol. XXVII (1903), p. 226. 
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§ 4. The Finite Integral of a Function which is Analytic in a Circular Ring. 

If $>{x) is any function which is analytic within a circular ring we may 
deal with the problem of solving the equation 

u(x + 1) — u(x) = <p(x) 

in a manner similar to that employed in the preceding section. "We suppose 
that ty{x) is analytic within the circular ring bounded by two circles each 
having the center a, the radii being r and R, where r + 1 < R. Take any two 
quantities f and R such that 

r <r, r + 1 < R < R, 

and form circles about a with radii r and R. Denote by r, r, R, R the circles 
having the radii r, r, R, R respectively. Take points A and B lying within 
the circular ring bounded by R and R and lying on a line through a parallel to 
the axis of iniaginaries, A being below and B above a. Let us define a region 
H in the following manner : Let the line AB cut r and R in the points L, I, m, M, 
these points being named in the order of increasing imaginary coefficients. 
From I and m draw tangents to r proceeding to the left and cutting R in the 
points T\ and T 2 respectively. The region H has the following boundary: 
The arc T 1 LMT 2 of the circle R, the tangent T 2 m, the right-hand semi-circle 
ml, the tangent IT X . Let »bea point in the region H and let X be the line 
through the points . . . . x — 1, x, x + 1, Form the function 

U ( X )= C ±M^_ (21) 

where the path of integration C extends from A to B, lies wholly in the ring 

bounded by r and R, cuts the line 7i in a point P between x — 1 and x, and has 

no other point in common with X. 

It is easy to show, by the method employed in the preceding section, that 

the function u(x), as thus defined, satisfies the relation u(x-\-l) — u(x) =q>(x). 

It is evidently analytic at every point of the region H. 

Now 

u(x — 1) = u(x) — q>(x — 1). 

By means of this equation the value of u(x) at a point x — 1 can be found 
when its value at x is known. In this way we may determine the function 
u(x) in that part of the circular ring between r and R in which it is not 
already defined. It is thus clear that u(x) is analytic throughout this circular 
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ring except (possibly) at the points congruent on the left to the singularities 
of <{>{x) in r. 

The right and left half -planes, as divided by the line AB, enter unsym- 
metrically into the above discussion. It is clear that we may make a similar 
discussion of the problem in which the two half-planes exchange roles. 

Thus we have the following two-fold theorem : 

V. If <p{x) is any function which is analytic within a circular ring 
bounded by two circles about a as center and having radii r and R where 
r + 1 < R, and if any other concentric circular ring is formed with radii 
r and R, where r < r, r + 1 < R < R, then there exists a function u(x) satis- 
fying the relation 

u(x + 1) — u{x) = q>(x) 

and having either of the following properties: 

1. It is analytic throughout the ring bounded by the circles f and R 
except at points congruent on the left to the singularities of $>{x) which lie in 
the circle r. 

2. It is analytic throughout the ring bounded by the circles r and R 
except at points congruent on the right to the singularities of <p(x) which lie 
in the circle r. 

II. Entire Functions Defined by Linear Homogeneous Equations with 

Rational Coefficients. 

Let us consider a single difference equation of the nth. order 

R (x)H(x + n) + R 1 (x)H(x + n — 1) + +R n (x)H(x) = $(x), 

where R (x), R x (x), . . . ., R n (x) are rational functions. Let a be the least 
integer such that x~ ia R i {x)/R G (x) is analytic at infinity for every value of 
i, i — 1,2, ....,n; and let R f be the constant term in the expansion of 
x~ ia R i (x) /R (x) in powers of x~ l . We assume that the constants i^ are such 
that the roots p 1} p 2 , . . . ., p w of the characteristic equation 

p» + R.f- 1 + i? 2 p»- 2 + .... + i?„_ lP + R n = 

are in absolute value different from each other and from zero. An equation 
having these properties, whether homogeneous or non-homogeneous, is here 
said to be of simple character. In what follows we confine attention for the 
most part to equations of simple character. 



176 Caemiohael: On the Theory of Linear Difference Equations. 

If the equation (whether of simple character or not) is homogeneous and 
the coefficients are constants it is easy to determine by elementary means a 
fundamental system of solutions, each function of which is entire. Guichard * 
has gone further and has shown that if the coefficients are constants arid the 
independent term $(x) is any entire function, then there exists a particular 
solution of the equation which is an entire function. In this case the general 
solution can be expressed in the form 

H(x) = E(x)-\- p, (x)E L (x) + .... + p n (x)E n (x), 

where E(x), E 1 (x), ....,E n (x) are entire functions and p 1 (x), ....,p n (x) 
are arbitrary periodic functions of period unity. 

The question naturally arises as to whether corresponding theorems are 
true for equations in which the coefficients are no longer constants, but are 
restricted to be rational functions. Such is the subject of investigation in this 
and the following section of the present paper. 

The theorem on page 133 of my dissertation, taken in connection with 
the remarks which follow, leads at once to a theorem which we shall need ; it 
may be stated thus : 

// 

P (x)H(x + n) + P 1 (x)H(x + n — 1) + . . . . +P n {x)H (x) = 

is a homogeneous difference equation of simple character in which the coeffi- 
cients P (x), P 1 (x), . . . ., P n (x) are polynomials in x, then 

1) There exists a fundamental system of solutions H x (x), ....,H n (x) 

each function of which is analytic throughout the finite plane except at the 

zeros of P n (x) and the points congruent to them on the left. If x approaches 

infinity in the positive direction along any line parallel to the real axis, we 

have 

lim H t (x) \T(x) \- a a j - x x- m > = 1, j = 1, , n, 

where a, a jf m^ are appropriately determined constants; 

2) There exists a fundamental system of solutions H 1 (x), . . . ., H n (x) 
each function of which is analytic throughout the finite plane except at the 
zeros of P (x — n) and the points congruent to them on the right. If x ap- 
proaches infinity in the negative direction along any line parallel to the real 
axis, we have 

lim H j (x) \~f (x) \ -*bf*x- n > = 1, j = l, , n, 



* Annates de VEcole Normale Supfrieure (3), 4 (1887), pp. 368-370. 

f Transactions of the American Mathematical Society, Vol. XII (1911), pp. 99-134. 
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where /?, b jf n i are appropriately determined constants and 

f(x) = (l—e 2 " ix )r(x), t=V^T. 

From, either of these fundamental systems of solutions of (1) we may 
obtain (at once and in an obvious way) a fundamental system which consists 
of entire functions; for this purpose, it is sufficient to multiply each of the 
given solutions by an appropriate periodic function of period 1 so that the 
zeros of the latter will balance off the poles of the former, thus giving a pro- 
duct function which has no singularity in the finite plane. 

Thus we see that all functions denned by linear homogeneous difference 
equations of simple character can be expressed in terms of entire functions 
and functions which are periodic of period 1. That the entire functions defined 
by an equation are sometimes not the most interesting ones which the equation 
may define is evidenced by the importance of the gamma function. It has 
been found of much more use than any entire function which satisfies the same 
difference equation as the gamma function itself. 

That non-homogeneous equations, even when the coefficients and the inde- 
pendent term are all polynomials, do not always possess a solution which is 
an entire function will be shown in part III by means of an example. 

III. Applications op the Pkeceding Results to Non-homogeneous 

Equations. 

The results of the preceding parts of the paper will now be applied to an 
investigation of properties of the solutions of non-homogeneous equations 
whose coefficients are polynomials in x but whose independent term is a mero- 
morphic function or an entire function of x. The results, other than the nega- 
tive one afforded by the following example, are stated in Theorems I and II. 

That non-homogeneous difference equations in which the coefficients and 
the independent term are polynomials do not always possess solutions which 
are entire functions is shown by the equation 

xg(x-\- l) — (x + l)g(x) = x, 

the general solution of which may be written in the form 

where p (x) is an arbitrary periodic function of x of period 1. It is obvious 
that g (x) has singularities in the finite plane for every choice of p (x) . But if 
p(x) is any periodic function which has in the finite plane no singularity other 
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than poles, then the corresponding g(x) is analytic throughout the finite plane 
except for poles. 

This example suggests the following theorem, which will now be proved : 

I. // 
P (x)H(x + n) + P 1 (x)H(x + n — 1) + . . . . + P n (x)H(x) .■=$(*) (1) 

is a difference equation of simple character in which the coefficients are poly- 
nomials in x, while ${x) is analytic throughout the finite, plane except for 
poles, then its general solution H(x) may be written in the form 

H{x) = M(x) + p l (x)E 1 (x) + ....+ p n (x)E n (x), 

where E x {x), . . . ., E n (x) are entire functions, M(x) is analytic throughout the 
finite plane except for poles, and p x {x), ...., p n {x) are arbitrary periodic 
functions of period 1. 

In section II we have shown that the corresponding homogeneous equation 
(that is, the equation obtained from (1) by writing zero for $>(x)) has a 
fundamental system of solutions each function of which is entire. We let 
Eyix), . . . ., E n (x) be the functions of this solution. It is now obvious that 
our theorem will have been demonstrated when we have proved the existence 
of a single function M(x) satisfying the non-homogeneous equation (1) and 
possessing the required character of being analytic throughout the finite plane 
except for poles. We shall now prove the existence of such a function. Our 
method is the classic one of variation of parameters. 

We seek a particular solution H(x) of (1) in the form 

H(x) = a 1 (x)E 1 (x) + a 2 (x)E 2 (x) + .... + a n (x)E n (x). (2) 

In order that H(x) shall be a solution of (1) it is sufficient that the functions 
a t (x) satisfy the set of equations 

E t {x + 1) Aa^) + .... + E\(x + 1) Aa n (a) = 0, 
E,(x + 2) A ai (*0 + .... + E n (x + 2) Aaja) = 0, 



> 



(3) 



E x {x + n — 1) A ai (x) + .... + E n (x + n — 1) Aa n (a?) = 0, 
# 1 (aj + ™)Aa 1 (a;) + + E n (% + w)Aa„(r) = ${x). 

This may readily be shown by the method usually employed in finding a solution 
by means of variation of parameters. It may easily be verified directly ; for, 
from (2) and the first equation in (3) one has 

H(x + 1) = a 1 (x)E 1 (x + 1) + . . . . + a n (x)E n (x + 1). 
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By means of this equation and the second one in (3) one obtains the result 
H(x + 2) = a 1 (x)E 1 (x + 2) + . . . . + a n (x)E n (x + 2). 

The process may evidently be continued till H{x-\-n) is found. Then a direct 
substitution in (1) will show that that equation is satisfied. 

The determinant of (3) is not identically zero; for this is the condition 
that E 1 (x), . . . ., E n (x) form a fundamental system of solutions of the 
homogeneous equation corresponding to (1). Hence (3) can be solved for 
Aaj(#), . . . ., Aa n (x). Thus we have 

Aa i (x) = M i (x), i = l, ,n, (4) 

where the functions M t (x), i — l,....,n, are obviously analytic throughout the 
finite plane except for poles. From Hurwitz' theorem, quoted above (p. 163), 
it follows then that there exist functions a^x), i = 1, . . . .,n, having the same 
property. If these functions a f (x) are put in (2) we have a function H(x) 
satisfying (1) and having the property of being analytic throughout the finite 
plane except for poles. We let H(x) be the function M(x) of the theorem. 
Thus the proof of the theorem is completed. 

We shall now turn our attention to a more restricted form of equation 
(1). We suppose that q>(x) is an entire function and that P„(a?) = l. We 
write 

P (x)H(x + ») + .... + P n -,{x)H{x + 1) + H(x) = <p{x). (1') 

Then according to the theorem quoted in part II above, there exists a funda- 
mental system of solutions, 

H,(x), H 2 (x), , H n (x), 

of the corresponding homogeneous equation each function of which is analytic 
throughout the finite plane ; and these functions verify the relations 

lim H f (x) \T(x)\- a p-*x~ m i = 1, j.= 1, , n, (5) 

where the limit is taken for x approaching infinity in the positive sense along 
any line parallel to the real axis. Here the constants p,-, j = 1, . . . ., n, are 
the roots of the associated algebraic equation and m ; , j = 1, . . . ., n, are prop- 
erly determined constants. 

We seek now a solution H^x) of (1') of the form 

H\x) = aMH.ix) + .... + « n (x)H n (x). 
23 
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Corresponding to equations (3) we have now the system 

H x (x + 1) A ttl («) + .... + H n (x + 1) Aa.(<*0 = 0, 

> 

H x (x + n — l)Aa 1 {x) + .... + H n (x + n — l)Aa n (x) = 0, 
^(a; + n)Aa x (x) + +#„(# + n)Aa n {x) = <p(x). 

The determinant of this system is D(x -f- 1) where 

#,(*) 5 2 (^) ....H n (x) 

H^x + 1) H 2 (x + 1) ....HJx + 1) 



(6) 



2)(aj) 



fl^aj + w — 1) jgr 8 (a? + w — 1) H„(a + w — 1) 



If We employ the notation lim ^(v) in the same sense as in (5), then in view 
of (5) and the definition of D(x) it is easy to prove that 

lim D(x)\F(x)T(x + l) 

Y{x + n — l)}~ a (p 1 p 2 p.) -•»-*-«— —«-=Z), (7) 

where 

1 1 1 

pr 1 p 2 _1 p^ 1 



D=l 



Pi P2 



P„" 



» + l 



P. 



whence ZHs a constant different from zero, since the quantities p lf p 2 , 
are all distinct. 

By Stirling's formula for the gamma function we have 

lim T(x + r) (x + r )- x - r +le x+r = V2n, 

where r is any integer. Hence we have readily 

lim T(xi- r)x- x e x x~ r +h=V2n. 

Taking this result in connection with (7) it is easy to see that there exist con- 
stants a, m, c, crf= 0, such that 

lim D(x)x nax a~ x a>~ m = c. (8) 

Now D(x) satisfies a difference equation of the first order which may be 
found in the following manner : In the determinant defining D (x) replace x by 
x + 1 and then substitute for H f (x + n) its value 



H t (x + n) =- -p-^- [P x (x)H i (x + n - 1) + . . . . +#,(#) ] 
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gotten by means of the homogeneous equation corresponding to (1). One thus 
obtains the result 

p(* + i, = (_i).£M. 

If we make the transformation 

D (x) = D (x) x na *a- x x- m 
we have readily the result that 

D(x + l) = \l+^{x)\D(x), lim D{x) = c, (9) 

where ^{x) is analytic at infinity and is expansible in the form 

>i(x) = V'x- 2 + 4>"'x~ s + , \%\>R, 

R being an appropriately chosen constant. That D(x) is uniquely defined by 
the fact that it satisfies relations (9) is easily shown. For, if D x (x) is any 
function satisfying the given difference equation then the most general solution 
of that equation is p(x)D 1 (x), where p{x) has the period 1. The existence 
and value of the limits lim D 1 (x) = c, lim p(x)D 1 (x) = c require that p (x) 
shall equal 1. Hence D(x) is uniquely defined by relations (9). 

In view of the general theorem in § 1 of my dissertation, p. 108, it is now 
clear that D (x) . is such that there exists a line parallel to the axis of imagin- 
aries to the right of which D (x) has no zero. Then, from the relation between 
D (x) and D (x) it follows that there exists a line I to the right of which D (x) 
has no zero. 

From these results we conclude that system (1') has a solution in the 
form 

Act i (x) = M i (x), i = l,,...,n, (10) 

where the functions M { (x) have no singularities to the right of the line I. A 
theorem of Hurwitz* then states that there exist functions a t (x) satisfying 
equations (10) and having the property that there exists a line X, parallel to 
the axis of imaginaries, to the right of which no one of these functions has a 
singularity. Hence, if we take such values of a { (x), the function H 1 (x) has 
no singularities to the left of the line "k. Then, making use of the fact that 
Wix) satisfies (1'), we conclude readily that H x (x) has no singularity in the 
finite plane. Hence the general solution of (1') may be written in the form 

H(x) = W(x) + p 1 (x)H 1 (x) + .... + p n (x)H n (x), 
where Px(x), . . . ., p n {x) are arbitrary periodic functions of period 1 and the 
functions 2J 1 (se), 2J 1 (#), . . - ., S n (x) have no singularities in the finite plane. 

*Acta Mathematica Vol. XX (1897), pp. 308-310. 
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In a similar manner it may be shown that if we modify (1) by taking 
P (x) instead of P n (x) equal to 1 and making $(x) entire as before, we obtain 
the same general result. 

Thus we have the following theorem : 

II. // 

P (x)H(x + n) + P 1 {x)H{x + n — 1) + . . . .+P n (x)H(x) = <p{x) 

is a difference equation of simple character in which $(x) is an entire function 

and the coefficients P i (x)., i = 0, 1, . . . ., n, are polynomials in x, and if either 

P (a;) = l or P n (x)s=l, then the general solution H(x) may be written in the 

form 

H(x) = E(x) + p 1 (x)E 1 (x) + ....+ p n (x)E n (x), 

where E(x), E- i {x), . . . ., E n {x) are entire functions and p±{x), . . . . , p n (x) 
are arbitrary periodic functions of period 1. 

Indiana University, April, 1912. 



